Minimal Multifractality in the Spectrum of a Quasiperiodic Hamiltonian 
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In many systems, the electronic energy spectrum is a continuous or singular continuous multifrac- 
tal set with a distribution of scaling exponents. Here, we show that for a quasiperiodic potential, the 
multifractal energy spectrum can have a minimal dispersion of the scaling exponents. This is made 
by tuning the ratio between self-energies in a tigth-binding Hamiltonian defined in a quasiperiodic 
Fibonacci chain. The tuning is calculated numerically from a trace map, and coincides with the place 
where the scaling exponents of the map, obtained from cycles with period six and two, are equal. 
The diffusion of electronic wave packets also reflects the minimal multifractal fractal nature of the 
spectrum, by reducing intermittency. The present result can help to simplify the task of studying 
several problems in quasiperiodic systems, since the effects of multiscaling are better isolated. 

PACS numbers: 



Quasiperiodic systems, which are neither periodic, nor 
disordered, have been the focus of an active research 
since the seventies 0, Q- Later on, they gained much 
more importance when the discovery of alloys with five- 
fold symmetry and long range order was made Q . These 
alloys were called quasicrystals, in recognition of their 
quasiperiodic atomic structure. The word quasiperiodic 
means that the structure can be represented as a sum 
of periodic functions, each with an incommensurate pe- 
riod with respect to the others. Still, the nature of the 
physical properties of quasicrystals is not well understood 
00 j^. Even in the theoretical side there is a lack of un- 
derstanding in how electrons propagate, specially in two 
and three dimensions 0. As is well known, a periodic po- 
tential satisfies the Bloch's theorem, which tells that the 
eigenstates of the Schroedinger equation are plane waves 
of delocalized nature, and the energy spectrum is contin- 
uous 8]. For disordered systems, like in the one dimen- 
sional (ID) Anderson model, all the states are localized 
corresponding to isolated eigenvalues ^ . In more dimen- 
sions, there is a mobilrh' edge which separates extended 
from localized states For most of the quasiperiodic 
systems in ID, the spectrum is neither continuous nor 
singular, instead a new type of spectrum, called singu- 
lar continuous is obtained [lOj . This kind of spectrum is 
similar to a Cantor set, and presents a multifractal na- 
ture. The corresponding eigenfunctions are called criti- 
cal, and also show self-similarity and multifractality. In 
two and three dimensions, the nature of the spectrum 
is not known, although there seems to be a kind of mo- 
bility edge 12j. However, even in ID, where large 
amount of work has been done, there are m any unsolved 
questions, like the nature of conductivity |14| or diffu- 
sivity ^3 ; the spectral statistics and the shape of many 
of the eigenfunctions For example, Machida 

et. al. tl8|| showed that for certain parameters of the 
quasiperiodic Harper equation, the distribution of level 
spacing follows an inverse power law. This tendency was 
explained by Geisel et. al. ^ ^ level clustering. 

More recently |20| , it has been argued that the clustering 



regimen of the distribution of level spacings in quasiperi- 
odic potentials was an artifact of an inappropriate way 
of making the unfolding procedure, and the fact that the 
distribution of gaps and bands were different was con- 
sidered as a signature of different scaling indices in the 
spectrum 0. 

Many of the difhculties that arise in the studies of 
quasiperiodicity, are due to the fact that all of the sys- 
tems that have been used up to now, present a mul- 
tifractal spectrum 0|[T9] 20] 21], i.e., the scaling with 
the system size is not given by only one scaling expo- 
nent (a), instead a distribution of exponents determined 
by a function called /(a) is usually found. This makes 
very difficult to relate the spectrum with the correspond- 
ing eigenfunctions, since they have a complicate shape 
due to multiscaling and intermittency. For example, the 
moments of a wave packet that is diffusing behave as 
< x'^{t) >~ t'^P^i) for a time i, where I3{q) is an exponent 
related with the spectral type. In principle, it seems 
that one needs to exclude any simple relation between 
spectral and diffusion exponents 15] 22]. However, as 
we show in this letter, the problem can be more trans- 
parent if the parameters of the quasiperiodic potential 
are chosen in such a way that the system presents a min- 
imal dispersion of scaling exponents. As a consequence, 
instead of dealing with a wide multifractal set, which 
makes all computations very demanding from a technical 
point of view, one considers a very narrow multifractal 
set, which is much closer to a pure monofractal. Thus, 
it is possible to improve the isolation of intermittency, 
without dealing with huge multiscaling. 

As starting point, we use the simplest model of a qua- 
sicrystal: the Fibonacci chain (FC). The FC is build as 
follows: consider two letters, A and B, and the substi- 
tution rules, A ^ B, and B AB. If one defines the 
first generation sequence a,s J-i = A and the second one 
as J-2 = BA, the subsequent chains are generated using 
the two previous rules, for instance, JF3 = ABA. Start- 
ing with an A, we construct the following sequences. A, 
B,AB, BAB, ABBAB, BAB ABB AB, and so on. Each 
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generation obtained by iteration of the rules is labeled 
with an index /. Is clear that the number of letters in 
each generation / is given by the Fibonacci numbers F{1) 
of generation I, which satisfy: F{1) = F{1 - 1) + F{1 ~ 2) 
with the initial conditions: F{0) = 1,F(1) = 1. A nat- 
ural model for a one dimensional quasicrystal is a tight- 
binding hamiltionian of the type, 



(1) 



where ipn is the wave-function at site n, t„ is the reso- 
nance integral between sites n and n -\- 1. V{n) is the 
atomic on-site potential and E arc the allowed energies. 
For the present purposes, t„ is set to 1 for all sites, and 
V{n) has two possible values, Va and Vb following the 
Fibonacci sequence. In fact, using a Fourier expansion, 
we can show that the Fibonacci potential is no more than 
a sum of well known Harper potentials. 



V{n) = V + 25V^V{s) cos (7rs0(2n + 1)) (2) 



where (/) is the golden mean </) = {^/b-\-l)/2,V = {Va/4>)~^ 
{Va/cI)^), and V{s) is the s harmonic of the Fourier series, 
V{s) = sin(7rs0)/s0. The strength of the quasiperiodicity 
is measured by the difference between site-energies SV — 
Va-Vb. 

The spectrum of this problem, i.e., the energies that 
satisfy eq.Q, are obtained by using the trace of the 
transfer matrices. These energies are those for which the 
trace of generation / (that we denote by 2xi{E)) of the 
transfer matrix satisfy — 1 < xi{E) < 1. In the case of a 
FC, xAE) can be obtained from a non-linear map given 
by Hi, 



xi+i{E) = 2xi{E)xi-i{E) - xi^2{E), 



(3) 



with the initial conditions X-i = 1, xq = (E + A)/2, xi — 
{E — X)/2, and A = \6V\ /2. Since the map contains three 
successive generations, one can consider the evolution of 
the trace as a trajectory in a 3D space. It has also been 
shown that the spectrum is a multifractal set |23| , and it 
was conjectured that the exponents describing the scal- 
ing of the spectrum near the upper band edge and band 
center were governed by cycles of periods 6 and 2 of the 
trace map. Using a linearized stability ana.lysis around 
these cycles, two scaling indices were found [23, one cor- 
responding to the period 6 cycle. 



ln{[l + 4(l + A2)2]i/V2(l + A2)}' 
In 



(4) 



and the other to the period 2 cycle, 

ln{8J-l+ [(8J-l)2-4]'/'}/2 




FIG. 1: Scaling exponents ae (solid line) and 0:2 (dotted line) 
as a function of SV. The arrow indicates the place where qs = 
Qf2. The numerical results are shown as circles. Inset: absolute 
value of the difference between tJmax — amin for diferent values 
of SV. The triangules, squares and circles are the numerical 
results for generations I = 9,12 and 15 respectively. The solid 
line is the absolute value of the difference between ae and a2 . 



with J = [3 -f (25 + 16A2)1/2]/8. in Ref. H , it was 
concluded that a2 was the minimum scaling exponent, 
related with the band- width of the upper band edges, 
oq was the maximal scaling exponent. Between them, 
a continuous distribution of exponents was found. That 
argument can be related with the shape of the wave- 
function, and in principle it is possible to conclude that 
wave-functions in the edges decay more rapidly than in 
the center of the spectrum Later on, it was proved 
that the 6 cycle does not always occur at the center or 
edges of the spectrum, and in fact, the question of which 
is the maximal contraction factor is an open question, 
since there are other cycles |2^. In spite of this, in fig. 1 
we show a plot of scaling exponents from cycles 6 and 2, 
where is clearly seen that not always 012 < ctQ , as was 
erroneously stated in ref. [2^ . 

Furthermore, there is a cross-over of the two scaling 
exponents. The exponent a2 is the minimum only from 
SV — to SV 10.5. The crossover occurs when a2 — 
uq . The special value of SV* where this happens can be 
calculated by making Eq. 10} equal to Eq. After 
many lengthy algebraic steps, it is possible to show that 
SV* is given by the irrational number. 



SV* = 2W13 



11 r- 



10.497929740830855719340651. 



a2 



h-ii 



(5) 



This value perfectly agrees with the plot of ag and a2. 
In the special value SV* , — a2 — d*j ~ 0.305. If 
as stated in ref. the 6 and 2 cycles were the ex- 

tremal scaling exponents, then amax = Qfmin at SV*, and 
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thus one should expect only one scaling exponent, corre- 
sponding to a pure monofractal. To test how far is this 
assumption from the real behavior of the spectrum, we 
have performed numerical calculations. First, the band 
widths of the spectrum were obtained by the iteration 
of the map Q, using quadruple precision (32 significant 
digits) for each value of 6V. The calculations were made 
for chains of generations Z = 9, 12 and 15. The max- 
imal and minimal scaling exponents have the following 
expression. 



\nAEo - \nF{l) 



\nAEo-\nF(l) 
lnA£;„,ax(0 



(6) 



where AEq is the initial band width, and AE^i^(l) 
(Aii^,„ax(0) the minimal (maximal) bandwidths of 
a / generation FC. The contribution of AEq tends to 
zero in an infinite system, but is needed in the present 
finite case. Figure 1 shows the numerical values of amax 
and amin for / = 15 compared with the theoretical val- 
ues. In spite that ag and a2 are not the maximal and 
minimal scaling exponents, they are close to the numer- 
ical values, and a minimal dispersion of the spectrum is 
seen at 6V* . The reason of the very close values of the 
real extremal exponents, amax and amin, compared with 
ae and 02, is because these cycles dominates the spec- 
trum when 6V — > 00. Notice that beyond SV — 12, there 
are no numerical points due to the increased difhculty 
in computing the band widths, since a lot of precision is 
required as 5V — > 00. 

It is also possible to calculate the multifractal distri- 
bution /(a) of scaling exponents. Let us denote the 
bandwidth of the i— th band of the spectrum by AEi. 
Each band contains the same measure or density of states 
1/F{1). Each band scales as: {AE,)°' = 1/F{1). A parti- 
tion function is defined by p7| . 



Tl{q,T) 



F{1) 

E 



1 



Fil) 



1 



AE, 



(7) 



This function has a limit when I goes to infinity, which is 
either zero or infinity, unless r and q are chosen in an ap- 
propriate way 27] such that Ti{q,T) = 1. This condition 
determines a function T{q). The fractal dimension for a 
set of points with the scaling a, is obtained by a Legendre 
transformation /(a) = —T{q) + aq, where a = dT{q)/dq. 
In figure 2, the evolution of /(a) versus a is shown for 
different values of 5V using generation 1 = 9. Notice how 
the width of the curves decrease and in fact is a minimum 
as SV SV*. The width of the curve f{a) = also de- 
fines the values of amin and amax • This can also bee seen 
in the inset of fig. 1, where we plot the absolute value 
of the difference Aa = |amax — amini for different gen- 
erations. Observe that /(a) has a finite width at SV*, 
and thus amin 7^ amax- This rules out the possibility of 
having a pure monofractal spectrum, which was some- 
thing to be desired. However, the numerical calculations 
suggests that SV* is the value where the spectrum has a 
minimal multifractal dispersion. 



a 0.4 




FIG. 2: Evolution of f{a) as SV is changed for a 
FC of generation / = 9. From right to left, SV = 
3, 4, 5, 6, 7, 8, 9, SV* , 11, 12. The special value SV is indicated 
with an arrow. As an example, Qmax and Ofmin are indicated 
for SV = 3. 



Since there is a close relationship between band scaling, 
localization and the dispersion of electronic wave func- 
tions [23, one can expected to observe the signature of 
minimal multifractality in the electron diffusion. A way 
of characterizing the propagation consists in the s tudy 
of the wavefunction moments as a function of time |l5| . 
which are defined as, 

F{1) 
n=l 

where j is the site position in the chain. ^{n,t) is the 
wavefunction at time t on site n, that evolves following 
the time dependent Schroedinger equation, in which E 
is replaced by the operator —id/dt . As was already 
stated, in quasiperiodic systems {x'^{t)) follows a power 
law behavior of the type (x^(t)) = T>{q)t'^^'-'^\ where I'(g) 
is a constant, known as the diffusion constant for the 
periodic case. In the present work, we have performed 
numerical simulations to obtain (3{q) by using chains of 
generation I = 15. As initial conditions, a delta function 
wave-packet is drop at time t = at site no. Differ- 
ent initial sites were used, but the results are almost hq 
independent. The numerical computations were made 
by solving the time dependent Schroedinger equation by 
two methods that gave similar results. The first was a di- 
rect numerical resolution of the differential Schroedinger 
equation using a discretization of time, and the second 
consists in a decomposition in cigcnfunctions. Since each 
eigenfunction evolves independently with time, the re- 
sulting wave function at any time is just a linear com- 
bination of eigenfunctions, where the coefficients of the 
combination are the projections of the initial condition 
into each eigenvector. For long times, the second method 
is much more accurate. Finally, the scaling exponents 
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FIG. 3: /3(g) as a function of the difference in self-energies 5V 
for q — 2,4:,6 and 8. A rigorous bound for /3(g) is shown with 
a grey line. Also, the scaling exponents ae and Q2 are shown 
with solid and dotted lines respectively. 



/3(g) were obtained by fitting a power law in a log-log 
plot of {x'^{t)) versus t. In all of the cases, the fit was 
made from t = 1000 to t = 50000 (all units have been set 
to one) in order to avoid initial transients. 

Figure 3 shows the values of /3(g) as a function of SV 
for the first eight moments. Observe that odd moments 
are not displayed since they are zero. The graph shows a 
very interesting feature which is a nearly collapse of /3(g) 
for g = 4, 6 and 8 at SV* . This suggests that /3(g) ^ /3* « 
0.16 for such values of g. A local maximum is obtained 



for g = 2. Since (3* is nearly g independent, and /3(g) 
is different from one, is clear that the wave functions 
show intermittency but almost no multiscaling, although 
this effect is mitigated since /3(2) has a lower value. In 
figure 3 we have also included the values of amin and 
ttmax as a comparison. Notice the crossover near the 
collapse of /3(g). The fact that /3(g) and a are represented 
in the same graph is not just a coincidence, since the 
scaling of the spectrum is related with wave-localization. 
In fact, the present work suggests that in a very rough 
approximation « /3*/2, which is similar to the result 
obtained from a analysis of the relationship between wave 
localization overlap and the spectrum |29|. Also, in figure 
3 we present a comparison with a strict mathematical 
bound that has recently been proved |30l | for /3(g) when 
SV > 8.0. The bound works as predicted, although is 
higher than the actual values of /3(g). 

In conclusion, we have shown a value of the Fibonacci 
potential for which there is a minimal dispersion in the 
scaling exponents of the multifractal spectrum. The 
value of the potential seems to de determined mainly by 
the six and two cycles. The behavior of wave packet 
diffusion refiects this fact by showing almost no multi- 
scaling. We hope that the present work will bring much 
more understanding into the properties of quasiperiodic 
Hamiltonians, since one can better isolate effects of mul- 
tiscaling. For example, it would be useful to study the 
evolution of the spectral statistics near this minimal con- 
dition. Since the trace map is similar for phonons or 
the case of quasiperiodic hopping integrals, the presented 
property can be extended to such cases, or different trace 
maps |3ll |. 
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